General conditions for quantum adiabatic evolution 
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Adiabaticity occurs when, during its evolution, a physical system remains in the instantaneous 
eigenstate of the hamiltonian. Unfortunately, existing results, such as the quantum adiabatic the- 
orem based on a slow down evolution (H(et), e — > 0), are insufficient to describe an evolution 
driven by the hamiltonian H(t) itself. Here we derive general criteria and exact bounds, for the 
state and its phase, ensuring an adiabatic evolution for any hamiltonian H(t). As a corollary we 
demonstrate that the commonly used condition of a slow hamiltonian variation rate, compared to 
the spectral gap, is indeed sufficient to ensure adiabaticity but only when the hamiltonian is real 
and non oscillating (for instance containing exponential or polynomial but no sinusoidal functions). 



PACS numbers: 03.65.Ca, 03.65.Ta, 03.65.Vf 

The "adiabatic" process, from the Greek a— (a-), not, 
cUa (dia), through, {ioLivw (bainen), to pass, was in- 
troduced by Carnot (in 1824) and W. J. M. Rankine 
(in 1858) in thermodynamics, then by Boltzmann (in 
1866) in classical mechanics [lj]. In 1928, Fritz London 
applied adiabatic process in chemical kinetics. Con- 
cerning the quantum physics, in 1911-1916 Paul Ehren- 
fest used adiabatic invariance in the development of the 
'Old Quantum Theory' and in 1928 Born and Fock 0] 
demonstrated the quantum adiabatic theorem. By def- 
inition, quantum adiabaticity occurs when, during its 
evolution driven by an hamiltonian H(t), a quantum 
state \^(t)) prepared in an eigenstate |n(0)) remains 
close to the instantaneous eigenstate \n(t)) (with a proper 
phase choice) as time t goes on. The basic concept 
of adiabaticity in quantum theory has been widely ap- 
plied in both theories and experiments. Applications 
range from energy level crossings, such as Landau- Zener 
transition, Born-Oppenheimer molecular coupling, col- 
lisional processes, quantum control or adiabatic quan- 
tum computation [J 0] • Unfortunately, even for the two- 
level system, no sufficient conditions are known to effi- 
ciently describe an adiabatic evolution driven by a gen- 
eral hamiltonian H(t) For instance, an example as 
simple as the Schwinger's hamiltonian (solved hereafter) 

6] Hit) = ^ ( . c ° s Lt sinee ~T), proves that neither 
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the "usual" adiabatic phase evolution f Q E n /K— i(n\n) 
nor the commonly used approximate adiabatic criterion 
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statement may look surprisiry 
presented in textbooks [3, [l 
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are sufficient (or necessary) to ensure adiabaticity. 

"tiGl butis 



It is indeed well known 
in NMR or in quantum optics (through Rabi oscillation) 
that resonant terms can lead to population transfer, i.e. 
to a non adiabatic behavior. This is linked to branch 
points, connecting the different eigenstates of the adi- 



abatic Hamiltonian, and explains for instance that non 
adiabatic behaviour exist when several successive transi- 
tions between pairs of levels occurs 
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condition (jT|)is not valid globally. 

It is therefore important to derive general conditions, 
for a system and its phase evolution, which ensure adi- 
abaticity. This is the goal of this article. As a corol- 
lary we will answer the still pending (even in the two- 
level case) question: why and when the standard con- 
dition ([T]), of a slow hamiltonian variation rate, com- 
pared to the frequency associated to the spectral gap 
E n I , is a sufficient adiabatic con- 
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dition. Indeed, we show that condition ([1]) is sufficient 
to ensure adiabaticity but only when the hamiltonian is 
real and non oscillating. 

Because almost all existing results, as the adiabatic 
criterion JT]) are based on the so called adiabatic limit of 
a slow down evolution, we shall first start by studying 
the standard results and by explaining why the standard 
adiabatic theorem can not help to solve the problem. 
Hopefully this part will also clarify the recent debate 
concerning the adiabatic phase and adiabatic criterion 

BmmmmmmiMmmmmmmis 

followin g th e (over-subtle) "rediscovery" by Marzlin and 
Sanders [3l| that condition fl} is not a sufficient one. We 
shall then derive exact bounds for adiabaticity. We then 
discussed their validity in a the general two-level case 
and their simplification in the case of a non oscillating 
hamiltonian. For clarity some lengthly calculations are 
reported in an appendix. 



STANDARD RESULTS 
Quantum adiabatic theorem 

The Born and Fock's quantum adiabatic theorem has 
been rigorously demonstrated, several times and by sev- 
eral different methods (see for instance 0, 0, 0, H3] and 
references therein), extended to the infinite dimensional 
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setting by Kato [33J, studied as a geometrical holonomy 
evolution by Berry (34[, extended to degenerate cases 
(without gap condition) [35| and to open quantum sys- 
tem [36j |. 

In the non degenerate (E m ^ E n ) case, the adiabatic 
theorem stipulates that: 

|¥ 6 (t)> - e -* /o*(^»/A-<<«-l*-> ) iTieft)) = 0(e) ► 0, (2) 

£—►0 

where evolution speed is controlled by e and the subscript 
stands for the H e (t) — H(et) evolution jjsBJ . Here the 
dot designates the time derivative and J * / = f Q f(t')dt'. 

To illustrate the limited practical utility of the theo- 
rem, let's suppose that an external laser field, with con- 
stant angular frequency Ri(t) — lu, is applied to a two- 
level system that we want to adiabatically drive by exper- 
imentally modifying two parameters: the coupling Rabi 
frequency (proportional to the square-root of the laser 
intensity) R^it) — f2(t), ana ^ ^ ne detuning of the laser 
from resonance i?3(t) = S(t). The hamiltonian is, in the 

rotating wave approximation: -| ( s7(t)e*"* °S(t)-a> ) = 
H(R 1 (t),R 2 (t),R 3 (t),R 4 (t)). Due to the R 4 (t) = tut 
term, slowing down the time would lead to ui(et) = (ojejt. 
When e — > 0, this would require reducing lj to zero which 
is experimentally impossible. Moreover, even in the static 
field (u> — 0) regime, the theorem applies but only if S 
and £1 can be slowed down simultaneously. The theorem 
says nothing about the adiabaticity if S(t) and f2(t) are 
varied independently with time. 

Although undoubtedly of great theoretical interest, as 
in the quantum adiabatic computation using interpolat- 
ing hamiltonian [37j], the theorem describes an evolution 
driven by H(et) with e — > and is obviously of no utility 
concerning the evolution driven by H (t) itself, as in this 
case e = 1 and cannot be reduced to zero. The theorem 
is then better formulated within the parameter domain 
than within the time domain 0, |38| : an evolution driven 
by H(R(i)) is adiabatic if the parameter path, between 
an initial R; n parameter value and a final one Rfi n , is 
followed infinitely slowly. 



Approximate adiabatic condition 

Contrary to the quantum adiabatic theorem, the ap- 
proximate adiabatic condition |T]) can be applied to H(t) 
itself. The origin of condition (fTJ) arises Q from the fact 
that the error term in Eq. ([2|) can be written [32| as 

O(e) = Yjm^n H e!!-e'1 + 0(e 2 ), wner e the linear e de- 
pendence is here only implicit and, deliberately but con- 
fusingly, hidden in |n e ). This has been the source of 
confusion [1 H EE EH, EI El when used with e = 1 
where \m e ) — \m). The confusion occurs because, even if 
derived without any proof by using e = 1 -t± 0, the crite- 
rion (JTJ) ensures an adiabatic evolution in almost all the 



known examples: Landau-Zener(-Stuckelberg), Rosen- 
Zener-Demkov, Nikitin, Zhu-Nakamura models or in the 
Rapid Adiabatic Passage or STImulated Raman Adia- 
batic Passage (STIRAP) processes, ... 0,0, H3- Impor- 
tant enough, as we shall see, all these examples use non- 
oscillating (exponential or polynomial) functions. There- 
fore, the simple idea of adiabaticity, given by the con- 
dition (fT]), of a small but finite variation rate of H(t) 
(compared to the spectral gap), is broadly used. Simi- 
larly, as extracted from equation (|2|) without any proof 
by using e = 1 0, an adiabatic phase evolution 



of J Q E n /h — i(n\h) is widely used [25j. However as 
mentioned in the introduction the Schwinger's example 
demonstrates that this "usual" adiabatic condition, as 
well as this "usual" adiabatic phase, is neither sufficiently 
nor necessarily to obtain an adiabatic evolution. 

To avoid any confusion the term 0(e 2 ) has to be eval- 
uated. This can be done for instance by giving an exact 
bound 32] on the adiabatic fidelity K^n)! such as j5(| 
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Similar bounds 0, HI, [4(| exists. They can be used to re- 
store the usual theorem [55] because = e^^ = e^ 
vanishes when e — > 0. However, they have severe limita- 
tions because, due to the integral term, they require a 
maximal evolution time T to provide an adiabatic evolu- 
tion when none is needed. This can be easily seen from 



the Schwinger's H(t) 
ample with u) = 1 s 



hiouj 



cost 



sin tie 
- 1 and 9 = 0.01. 
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GENERAL BOUNDS 

In order to derive a more useful bound than ©, let's 
study the evolution of \^(t)) driven by a general N- 
level hamiltonian Hit). For the corresponding eigen- 
values E m (t) of H(t), the eigenvectors e l9m ^ \m(t)), 
m = 1,- • ■ , N form a so called adiabatic basis, where 
9 m {t) are arbitrary phases to be chosen conveniently 
later. To study the adiabatic evolution we assume that 
\^>(t = 0)) = \n(0)). The Schrodinger equation for 
l*(*)) = EH=i U mn (t)e i8 ^\m(t)), i.e. with U(0) = I, 
leads to the time-evolution equation: ihU = H'U , where 



H' mk = (E m + M m )5 mk - th(m\k)e l( - e 



(4) 



As usual, we identify the operators and their matrices 
in the standard (also called natural or canonical) basis 
|m st ), m = 1, • • • , N. Thus, H' = P~ l HP - ihP^P 
where the columns of P are the eigenvectors e l8m ^\m): 



\P\k s 



JO h (t), 



m st \k) of if. 
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The evolution is adiabatic if and only if the fidelity 
\U nn (t)\ is close to unity (or |||*)(*| - \n)(n\\\ < 1). In 
order to also study the phase evolution of the state we 
compare the matrix U to another time evolution matrix 
U' which can be more easily evaluated. 

Let's define U 1 by J7'(0) = 1 and 



and condition J8l) becomes: 



(m\n) 



iMJ' = p'- x h'p'u', 



(5) 



where P' is an auxiliary matrix to be chosen conveniently. 
Then, the important equality comparing two operators 

U(t)-U'(t) = (P'(t)-l)U'(t)-U(t)(P'(0)-l)- 

U(t) f U-^t^P'it'p'it'W (6) 
Jo 

can be established by multiplying it by U^ 1 and then 
taking the time derivative. 

Several choices are possible but, for simplicity we 
choose P' to have P'^H'P' as an eigenvalue value de- 
composition of H 1 . In this case U' is diagonal: U' nn = 
e - J / )/ h dt where E' n is the eigenvalue of the n th 
eigenvector \n') = P'|n st ) of H' . P' is unitary, so 
||P'|| = \\U\\ = \\U'\\ = 1. We then apply Eq. © on 
|n st ) and take the norm on both sides to have 



At)) 



<|||n'(0))-|n st )| 
\n')\\ (7) 



^(E n -E m )/h-i(n\h)+i(m\rh)-f t 

(12) 

This condition, first derived in [201 ] . generalizes condi- 
tion (JT|) when H is not real57J. However, as |T]), it 
is an insufficient adiabatic criterion for two reasons: it 
arises from a perturbative approach, and it neglects the 
condition ([5]), which is important for oscillating H. In- 
deed, it is only when the hamiltonian matrix elements 
are non oscillating functions - in the approximate sens 
of none of their sum, product, division or combination 
has a large number of monotonic changes - that the con- 
dition © can be neglected. More precisely, in the gen- 
eral case when all the P' nn — (m st \n') are real (or with 
a time independent phase argument) and monotonic, an 
important simplification occurs because J Q \P' nn (t')\dt' — 

\pL(t) - plm\ < \p; nn (t) - 1| + |p; rm (o) - 1|. in 

this case, we see, by using the 1-norm [56] , that the 
derivative condition © essentially reduces to the sole ([5]) 
condition. Similarly P' mn piecewise functions with finite 
number (A/ — 1) of monoticity changes [58j would lead to 
a E m \ p Ln - 1| < 1/M type of condition. 



Multi levels system 



Eq. ([7]) gives a bound as well as the correct phase evo- 
lution for adiabatic evolution. Sufficient adiabatic condi- 
tions are: 



We use here an exact perturbation theory [4l| to cal- 
culate P'\n st ). We write H' = H + V where V is a 
perturbation. For simplicity, i.e. in order to isolate the 
subspace, we renumber the states to have n = 1 and, 



\n'(t)) - \n s 



< 1 



(8) 



dt' = \\\n'(t'))\\dt' < 1 (9) H' = H„ + V; H Q = 



\(P'(t)-l)\n si )\ 
' P'(t')\n st ) 



To tight these bounds we choose the phase of \n!) to 
be such that (n st \n') > 0. The adiabatic fidelity is 
bound by the inequality 2(1 - \(^(t)\n{t))\) < |||#(t)) - 

ft E 'n 

e~ l Jo ~\n(t))\\ , which should now be tighten as much 
as possible by choosing the 8 m (t) phases. 

Links and differences, of Eq. ([7]) with the usual the- 
orem given by Eq. ((2]) and of Eq. ([8]) with the usual 
condition given by Eq. {T]), can be inferred by applying 
standard perturbation theory to Eq. (j4|): 



using the 1 + (N 
(see Eq. g])) 



1) block matrix notation, we choose 



H 'nn 

H',-h8' 



V = '1 ( o n't 
2 V« o 



We then apply techniques, detailed in the appendix, to 
endup with the following simple conditions: 



Ifi'l 



dt 



(S 1 



HOlllfi'll « 1 (13) 
d 



dt 



-n' 



< 1 (14) 



which are together sufficient adiabatic conditions because 
they imply Eqs. © and ©. 

Eq. (TTJJ is here to prevent the use of oscillating hamil- 
tonian. Indeed, as discussed previously, if H is real and 



E n w E n - ih(n\h) + h9 n + ^ Rl —ft, (10) "non-oscillating", meaning that Q' and (5 



H' — H' 1 



ik 



(11) 



Using the equality 9 m = 0„ + arg(— i(m\h}) for all m ^ n, 
creates (to this second order approximation) reals P mn , 



(piecewise) real monotonic functions, the condition (TTj 
essentially reduces to condition (jT3j) . 

Eq. (fl"3"|) itself can be seen as a generalization of the Eq. 
(TT2~j) which itself generalizes the standard condition fl}. 
Indeed, if we add to the condition (TIB")) the, fortunately 
common condition of negligible coupling within the space 
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orthogonal to |n), i.e. negligible [59j 5' off diagonal terms 
to have (6'~ 1 ) mm w (<5'mm) _1 , we can recover Eq. (TH?)) 
by choosing 6? m = 0„ + arg(— i(m|n)) (i.e. f2' real). 

Finally, the appendix indicates that, for a strongly 
non-oscillating (very few monotonicity changes) real 
hamiltonian H the sole usual condition ([1]) , which is then 
condition (|13p , is sufficient to ensure an adiabatic behav- 
ior. 



Two level system 

Let's now illustrate the results in the two-level (N = 2) 
framework. We write, by removing the average di- 
agonal energy, the general hamiltonian in the (spin- 
magnetic interaction H = —jB.^a) form: H(t) = 

huj (t) t cos0(i) sine(t)e" iv(t) 
2 Vsin^tyM - C os0(t) 

6i + arg(— i(2|l)), the appendix shows that Eq. 
Eq. (|9|) are equivalent to: 



Using —6\ 



| up sin 9 



ub cos 9 



w o - ft ar s(^ sin _ 



dt' 



= H1 

|<5'| 

d \n'(t')\ 



dt' S'(t') 



and 

< f (15) 

< 1 (16) 



In order to check their validity, or similarly the one 
of conditions (fl3|) and (fT4l) . we first use the simple ex- 
ample due to Schwinger [6j, where all the parameters 
uj , 9,(p = lo are real and time independent. In this case 
the condition ([TBI ) vanishes and U(t) = er l $a H l h = 



(COS -f Ijy- sin -f-) 



(COS — 2 ~f"*Ti — Sm 



where £Ir = 



-\/|ri'| 2 + <5' 2 is the generalized Rabi frequency. The adia- 
batic evolution (negligible off-diagonal terms in U) is en- 



sured by the condition -U-J- 



|n' | _ 



< 1 



n R ^J{u -u cos 0) 2 +u 2 sin 2 6 

which is indeed equivalent to our condition (115j) : = 

-, — l" slnfl <^ i Furthermore, our equation d). includ- 

ing its phase / -Ej/fi. = flnt/2, correctly describes an 
adiabatic evolution. 

On the contrary, using this analytical example (by 
looking at the resonant lo ~ lo or small 9 cases for in- 
stance) it is straightforward to demonstrates that Eq. 
©■ j^j = ^iL'o^ ^ 1j as wen as the "usual" adiabatic 
phase evolution uj t/2 = f Q Ei/h - i(l|l) (see Eq. ((2J), 
are not correlated with an adiabatic evolution. 

We now add to our study the condition (fl~6|) by the 
use of the real cycling hamiltonian H = 5 ( $ £5 ) where 
8{t) = acos(nrf) and a,w, f2 are positive constants veri- 
fying, for simplicity, weak-coupling (il <C a) and large 
amplitude (a > to). For t € [0, Ti = tt/to], ^ = 



change, so our second condition (Tl6[) reduces to condi- 
tion (|15p . The non- adiabatic transition probability pi 
(so called single-passage or one-way transition), is given 

by the Landau-Zener's formula: k e~f ™ [42! ] and the 
adiabatic limit pi — > is covered by the condition (|15p : 
= ffr < 1- After M (even) multiple 
passage, for t = MTi, the non-adiabatic transition prob- 
ability becomes pu ~ Pi s ™ os i' 1 ^' and depends of a relative 
(Stiickelberg) phase ~ — of the wavefunctionji^. For 
~ 7r/2[7r], j»m can be M times higher than p\ leading 
to a full non adiabaticity pm ~ 1 even if p\ <C 1 . This 
illustrates why, in such an oscillating case, condition (|15|) 
(fp^ ^ 1 ) is not sufficient and the extra condition (flU)) 
(7^ "C 1/M) is needed to ensure an adiabatic evolution. 
This example shows that, with an oscillating hamilto- 
nian, even if a single passage is quasi-adiabatic construc- 
tive interferences might accumulate the small non adia- 
batic amplitude to result, after multiple passages, in a full 
non-adiabatic transitioned]. This is very similar to the 
case of single crossing but with several levels Til 43j. or 
to multilevel system [19( , leading, using stationary phase 
(saddle-point) theorem or steepest descent WKB type 
of methods, to sums or products of dephased Landau- 
Dykhne-Davis-Pechukas's formulas corresponding to sev- 
eral successive transitions between pairs of levels [171 1 1 81 ) . 
Finally, this shows that the standard condition ([1]) breaks 
down, not only when resonant terms are present, as some- 
times believed [Ly, |2jJ, |22|, |44| , but more generally when 
oscillating terms are presents. 



\o\ 



(<5 2 +n 2 ) 3 / 2 



is real and with a single monotonicity 



CONCLUSION 

By simply diagonalizing the hamiltonian H' (hamilto- 
nian in the adiabatic basis), we have derived simple con- 
ditions, Eqs.([T3]) and (HU and exact bounds (Eq. Q) for 
the state and its phase, ensuring an adiabatic evolution. 
The usual (or standard) condition |T]) is found to be a suf- 
ficient adiabatic condition but only for a real and "non- 
oscillating" hamiltonian evolution. This explains why all 
the previously cited examples (Landau-Zener, STIRAP, 
...) deal with the (real) interaction representation or the 
dressed state basis, where to = 0, and use non oscillating 
functions such as exponential or polynomial ones. 

Condition (Ti"4|) prevents oscillation [6l| but unfortu- 
nately with no distinction between case with constructive 
crossings or case with destructive (Stiickelberg) interfer- 
ences. However, the generic most common case concerns 
a "complex enough" system with small total probabil- 
ity when the single crossing probability is small [45| . i.e. 
where the sole Eq. (fl~3|) . or Eq. ([I]) for real hamiltonian, 
is sufficient to ensure an adiabatic evolution. 

This result simply highlight the fact that the standard 
mathematical technique (so called asymptotic analysis) 
to study the adiabaticity consists in extracting, form the 
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global solution of the Schroedinger equation, a set of lo- 
cal solutions which individually covers a region (let say 
between time and T) , with a controlled behavior of the 
coefficients in the equation. This means that the crite- 
rion JT]) is local and that in order to study the adiabatic 
behavior of a given hamiltonian, one should cut its evo- 
lution in part where we could apply safely the criterion 
([1]), namely in part with single branching point or with 
single crossing between pairs of levels. Globally we shall 
add each local non-adiabatic amplitu de to g et the global 
non-adiabatic amplitude 



m 12, i 



9|, |43(. We would 
stress that all this should be very well known, but seems 
to be forgot by many physicist if we refer to recent pub- 
lished articles. Our article, demonstrate in a simple way 
that using non-oscillating function the number of local 
solution is obviously finite and so the added probability 
remains small if the criterion |T]) is globally fulfilled. 

Finally, the adiabatic evolution is strongly related to 
the (semi-)classical limit h — > of quantum mechanics 
[55] [li], to the WKB approximation HH, to the Mini- 
mal work principle [13] , to the quasistatic thermodynam- 
ical process (48|, and to perturbation theory. Therefore, 
we hope that this work and the given examples can en- 
able the development of significant techniques, or provide 
novel insights into these important systems. 

Thanks to Sabine Jansen to have pointed out to me 
Born and Fock's consideration concerning monotonicity. 



APPENDIX 
Multi levels model 

We demonstrate here that conditions (jT3|) and (| 14[) im- 
ply the conditions ((8]) and (J9j) . The techniques are similar 
to one used in Davis-Kahan sin 9 theorem or Weyl-Bauer- 



Fike's types of perturbative bounds [43, [5fJ, |5f|. The 
starting point is the exact Brillouin-Wigner perturbation 
theory that we demonstrate here for completeness 41 1 . 

We define the projector Q n = 1 — \n st )(n st \, which 
in matrix notation is Q n = ( [j 5 ) , and the eigenvector 
|n') oc | to') of H' = Ho + V with a simple normalization 
(n st |n') = f. Hq = H' — V commutes with Q n so^ 
(E' n - H )Q n \n') = Q n (E' n -H' + V)\n>) = Q n V|n'). 
When multiply by (E' n — Ho)^ 1 this directly lead, using 
Q n \n') = |n') — \n st ), to the Brillouin-Wigner equation: 
|n') = (1 - {E' n - Ho^QnV)- 1 ^). Using the matrix 
notation, and the blockwise inversion, this becomes: 



n' = 



(17) 



where HA' = E' n - H' nn = (n st \H + V\n'} 
(n st |U|n') satisfies 



_ TT I _ 

nn 



S'-iA' = ?— — (l + J'^A') 



The idea is now to use the smallness of 6' 1 f2' (see 
Eq. (Tl3]0 to evaluate |n') = |n')/^(n'|n') and its time 
derivative, i.e. to study Eqs. ([8]) and ([9]). 

We first take the norm of Eq. (fT8|) and use ||(1 + 
§>-i A >yi\\ < J2Z ll^-^'H* = (l-ll^'-iA'll)- 1 to see 
that: if H^'" 1 1| ||f2'|| < f then H^A'H < f (see also Eq. 
([15])). Therefore Eq. (JT7J) shows that Eq. © is implied 
by (it is in fact equivalent to) Eq. (T3J: ll^ -1 !! ||fi'|| <C 1. 



When ll^'- 1 
derivative of 



in'H < 1, i.e. 



the time 



to') = |n')/ v /(n'|n') shows that |||n')|| » 
|||n')||. Time derivative of the equations ([T7|) and ([T8|) 
can then be used to study Eq. ([9|). Indeed, using 
<5' _1 = — <5' 1 <5'<5' 1 and (again) useful estimations on the 
smallness (see Eq. (|18jl) of the norm of <5' _1 A' and its 
time derivative, finally leads to the fact that condition 
(together with ([T5)0 implies the condition ©. 



Two levels model 

We derive here, in a simpler way, the conditions 
([15]) and (H]). The eigenvectors e Wl |1), e if>2 12) of the 
hamiltonian H , corresponding respectively to the eigen- 
values fkuo/2 and — Tkjo/2, are given by the columns 

V*f cosfe i£, i -e- 4 * sinfe i£, 2 ' 



of P = 



and H' 



2 \((psm9 -i9)e l( - 9l ~ 9 ^ -u + <p cos 9 + 29 2 

With 9 2 = -9i, H' = | ( n 'g, V The N = 2 case is 
a very special one because it is always possible to choose 
H', and then P' real with 9 2 = 6 X + arg(-i(2|i)). Using 

the obvious notations H' = ^ ( cose ', sin9 ', ) , i.e. P' = 

2 V sin 9 — cos J ' 

Our conditions (JHKS]) then read Q l < 1 



cos — sin - 
sin ^- cos % 
■t i 



(18) 



and J Q \9 \ -C I which leads to the general conditions of 
adiabatic evolution: Eqs. ([T5]) and <\1Q\i . 



Non oscillating case 

We assume here that the usual condition ([1]) is fulfilled 
for a strongly non-oscillating real hamiltonian (i.e. with 
monotonies the P m k — (m st \k) = tp m k functions) and we 
give here a clue that the evolution is indeed adiabatic. 

By using a proof by contradiction, we assume that the 
evolution is not adiabatic. Thus condition p3[) is not ful- 
filled so non negligible S' off diagonal terms exists to mod- 
ify substantively the 6' eigenvalues. The (Weyl-)Bauer- 
Fike's theorem ([T9]) applied to 6' , implies that one of the 
off diagonal elements (~ ip m k) of 5' should then be bigger 
than the diagonal ones (the gap AE n ). But Eq. ([20]) in- 
dicates that a time T ~ l/||fi'|| is needed to have an non 
adiabatic evolution. Thus condition ([T]) , which is roughly 
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||f2'|| <C AE n , would implies that ip m k ~ <Pm,kT 1 
which contradicts <p m k = (m st \k) < 1. 

(Weyl-)Bauer-Fike's theorem 

Let H' d be the diagonal part of H' = H' d + H' non diag . 
Multiplying K ondias \n') = (E' n - H' d )\n') by (E' n - 
•H'd) -1 and taking norm on both sides leads to the (Weyl- 
Lidskii)-Bauer-Fike's theorem (applied to H'): 

min\E' n -H' mm \ < \\H^ on diag || (19) 



Universal optimal bound 

Using the following choice 9 m — J * i(m\rh) — f Q E m /H 
of a geometrical phase (Berry Phase for cyclic evolution) 
plus a dynamical phase simplifies the H' matrix elements 
(see Eq. gj). Using -^jM < |^f-| and the norm 

[56] equality y/1 - \U* n \ = ^/E m# J^ n | 2 = ll^nll, 

when integrating the (Schrodinger) equation (ihU nn — 
Q'U, n ) leads to the (quantum Zeno's type of) adiabatic 
condition: 

1 - \U nn (t)\ < 1 - cosdlO'Hi/2) < ^J^t 2 . (20) 

This optimal bound is reached by the Schwinger system 
for §' = 0. 
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